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ABSTRACT 

A geostatistical method for gradually deforming an initial distribution of 
objects, of geologic type for example, from measurements or observations, so as to 
best adapt it to imposed physical constraints of, for example, a hydrodynamic type 
having applications of geostatistical modelling of heterogeneous reservoirs of 
various objects: fracture, channels, vesicles, etc., for example.[-] The objects are 
distributed in a zone of a heterogeneous medium according to a Poisson point 
process in form of figurative points with a point density X{x) that varies according to 
their position (x) in the zone, a realization of a uniform random vector according to 
which the position of each object is defined while respecting density X(x) is formed, 
and the uniform random vector is gradually modified according to a gradual 
deformation process so as to obtain gradual migration of each object until a final 
realization best adjusted to parameters relative to the structure of the medium, 
such as hydrodynamic parameters, is obtained. 



BACKGROUND OF THE INVENTION 
Field of the Invention 

[0001] The object of the present invention is a method for gradually 
deforming an initial distribution of objects of geologic nature, formed by simulation 
of a stochastic model of object type, from measurements or observations, so as to 
best adapt it to imposed physical constraints of, for example, a hydrodynamic type. 

Description of the Prior Art 
[0002] French Patent 2,780,798 filed by the assignee describes a method for 

gradually deforming a stochastic model (of Gaussian type or similar) of a 

heterogeneous medium such as an underground zone, constrained by a series of 

parameters relative to the structure of the medium. This method comprises drawing 

a number p of independent realizations (or representations) of the model or of at 

least part of the selected medium model from all of the possible realizations, and 

one or more iterative stages of gradual deformation of the model by carrying out 

one o r m ore s uccessive I inear c ombinations o f p i ndependent i nitial r ealizations, 

then composite realizations successively obtained possibly with new d raws, etc., 

the coefficients of this combination being such that the sum of their squares is 1 . 

[0003] French Patent 2,795,841 filed by the assignee describes another 
method for g radually deforming the representations o r realizations, g enerated by 
sequential simulation, of a non necessarily Gaussian stochastic model of a physical 
quantity z in a gridded heterogeneous medium in order to adjust them to a series of 
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data relative to the structure or to the state of the medium, which are collected by 
prior measurements and observations. The method essentially comprises applying 
an algorithm allowing gradual deformation of a stochastic model to a Gaussian 
vector with N mutually independent variables, which is connected to a uniform 
vector with N mutually independent uniform variables by the Gaussian distribution 
function so as to define realizations of the uniform vector, and using these 
realizations to generate representations of this physical quantity z, which are 
adjusted in relation to the data. 

[0004] The above methods are applicable to gridded models (pixel type 
models) suited for modelling continuous quantity fields and they are therefore ill- 
suited for modelling of zones crossed by fracture networks or channel systems for 
example. 

[0005] Models based on objects are spatial arrangements of a population of 
geometrically defined objects. Basically, an object type model is a Boolean model 
that can be defined as a combination of objects identical by nature with a random 
spatial distribution. Boolean models (of object type) are of great interest for the 
geometric description of heterogeneous media such as meandrous deposit 
systems, fracture networks, porous media on the grain size scale, vesicle media, 
etc. Geologic objects are defined by their shape and size. Their location in the field 
is defined by taking account of their interactions : attraction-repulsion, clustering 
tendency, etc. 
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[0006] Unlike pixel type models, the models based on objects can provide for 
example realistic geologic representations of an underground reservoir at an early 
stage where the data obtained by in-situ measurement are still rare. 

[0007] The prior art in the sphere of object type models is notably described 
in the following publications : 

Matheron, G., 1967,"Element Pour une Theorie des Milieux Poreux", Masson, 
Paris ; 

Matheron, G., 1975,"Random sets and Integral Geometry", Wiley, New York ; 

Serra, J., 1982, "Image Analysis and Mathematical Geology". Vol. I, Academic 
Press, London ; 

Stoyan, D.S. et al., 1995, "Stochastic Geometry and its Applications", 2"^ Edition, 
Wiley, Chichester ; 

Lantuejoul, C, 1997, "Iterative Algorithms for Conditional Simulations , in Baafi and 
others, eds."; Geostatistics Wollongong 96, Vol.1, Kluwer Acad. Publ., Dordrecht, 
The Netherlands, p.27-40. 

[0008] The position of the objects in an object type model is distributed 
according to the Poisson point process. The shape and the size of the objects are 
independent of their positions. This model can be generalized by a combination of 
objects of different nature and/or using a non-stationary density Poisson point 
process. 

[0009] Although Boolean models have been widely studied in the literature, 
there is no coherent and efficient method for constraining these models to the 
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physical data, notably hydrodynamic data, which is however a major challenge for 
their application to reservoir engineering. The methods allowing gradual 
deformation of pixel type stochastic reservoir models such as those described for 
example in the two aforementioned patents cannot be directly used for Boolean 
models. Constraining Boolean models to hydrodynamic data for example requires 
development of coherent algorithms for deformation and displacement of the 
objects. 

SUMMARY OF THE INVENTION 

[0010] The method according to the invention allows generalization of the 
gradual deformation technique described in the two aforementioned patents to 
stationary or non-stationary Boolean models, with or without geometric constraint 
as regards the well. The method appears to be particularly useful notably for 
reservoir engineers anxious to coherently and efficiently adjust object type reservoir 
models. 

[0011] The method according to the invention allows gradual deformation of 
an initial realization defining the distribution of a set of objects in a zone of a 
heterogeneous medium such as a geologic structure, formed by simulation of an 
object type stochastic model, the objects being distributed in the zone according to 
a Poisson point process in form of figurative points with a point density X{x) that 
varies according to their position (x) in the zone. It essentially comprises the 
following steps: 
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generating a realization of a uniform random vector according to which the position 
of each object is defined while respecting density 7i(x); and 

gradually modifying the uniform random vector according to a gradual deformation 
procedure, so as to obtain gradual migration of each object and consequently 
gradual change in the distribution of the objects in the zone, until a final realization 
best adjusted to parameters relative to the structure of the medium, such as 
hydrodynamic parameters, is obtained. 

[0012] It is possible to limit the migration of the figurative points in a 
subdomain of the zone (a well crossing the zone for example) by imposing a zero 
point density in the complementary part of the subdomain. 

[0013] According to an implementation mode, one gradually changes from a 
realization containing a first set of Ni points to a realization containing a second set 
of N2 points by constructing a chain N(t) of Poisson numbers between the two 
numbers Ni and N2, using the gradual deformation procedure. 

[0014] It is also possible to gradually deform the size, the shape and the 
orientation of an object during its migration by using the gradual deformation 
procedure. 

[0015] In cases where it contains errors, it is also possible to gradually adjust 
the point density X{x) by using the gradual deformation procedure. 

[0016] The final realization obtained by means of the method is a realistic 
representation of the configuration of the objects in the heterogeneous medium. 
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[0017] The method according to the invention finds applications in the sphere 
of underground zones modelling intended to generate representations showing how 
a certain physical quantity is distributed in a subsoil zone (permeability for example) 
and which are best compatible with observed or measured data : geologic data, 
seismic records, measurements obtained in wells, notably measurements of the 
variations with time of the pressure and of the flow rate of fluids from a reservoir, 
etc. 



BRIEF DESCRIPTION OF THE DRAWINGS 

[0018] Other features and advantages of the method according to the 
invention will be clear from reading the detailed description hereafter, with 
reference to the accompanying drawings wherein : 

[0019] [-] Figs.lA, 1B respectively show a density function of a non- 
stationary Poisson point process and a realization of the non-stationary Poisson 
point process generated by the sequential method ; 

[0020] [-] Figs. 2A to 2D show various examples of gradual migration 
trajectories between two points of a non-stationary Poisson point process ; 

[0021] [-] Figs.SA to 3H show various successive stages of a realization of 
the non-stationary Poisson point process during the gradual migration of the 
points ; 

[0022] [-] Fig. 4 shows the possible domains of migration of a disk in the 
case of three conditioning points ; 
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[0023] H Fig. 5 shows, by way of example, a complete chain of successive 
realizations of a Boolean simulation of elliptical objects where deformation 
parameter t ranges between -n and n by intervals At=0.l7i ; and 

[0024] [-] Fig. 6 shows, by way of example, an incomplete chain of 
successive realizations of a Boolean simulation of elliptical objects where 
deformation parameter t ranges between 0 and 0.2n by intervals At=0.0l7c. 

DETAILED DESCRIPTION 

General points 

[0025] The geologic objects to which the method applies are for example 
more or less large fractures inside a reservoir zone, or channels. The method can 
also apply to granular or vesicle type structures of much smaller size. All these 
objects are difficult to model by means of pixel type models. 

[0026] The progressive migration operations described hereafter apply to an 
initial model where the positions of the objects are represented by point 
configurations (referred to as point procedures). The distribution of these points 
varies in density according to their positions in the modelled zone. This distribution 
is based on various data known from measurements or 
observations : geomechanical measurements obtained in wells for example, 
seismic data obtained from prior seismic operations. 
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[0027] Starting from this initial model, the method allows gradual deformation 
of the initial distribution according to a series of rules so that, in the final 
distribution, the model is optimized in order to better correspond to constraints 
imposed on one or more physical parameters such as, for example, a distribution of 
permeability values. The migration rules are such that a great number of different 
points of the model can be globally displaced from a limited number of control 
parameters. 

[0028] An objective function measuring the difference between the physical 
parameters from the real heterogeneous medium and those simulated on a 
realization of the stochastic model is generally defined. The value of the objective 
function therefore depends on control parameters of the stochastic model. The 
values of these control parameters are obtained by minimizing the objective 
function. 

Summary of the Poisson point process 

[0029] The Poisson point process is a countable random set of points 
distributed in any space This set of points has the following characteristics : 

[0030] Let D be a domain of 9?". If the volume of D, denoted by I D I, is finite, 
then the numbers of points lying in D, denoted by N(D), follows a Poisson's law of 
parameter X |D I. Let : 

P[N(D) = «>-^l^l ^^1^ V« > 0 (1) 
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where X is called the point process density ; it measures the average number of 
points lying in a unitary volume domain of 9?". 

[0031] Let Di, D2,..., Dk be domains of 9?" unconnected two by two, then the 
numbers of points lying in these domains N(Di), N(D2),..., N(Dk) are mutually 
independent random variables. 

[0032] Conditionally on N(D)=np, these np points are independent and evenly 
distributed in D. 

[0033] The d escription h ereafter p resents t he P oisson point p rocess i n the 
finite domain D. 

Migration of a stationarv Poisson point process 

[0034] Consider the problem of the migration of a realization of a stationary 
Poisson point process in D rectangular. In order to simplify the presentation, D is 
assumed to be a unitary hypercube [0,1]" with n dimensions. Let Xi and X2 be two 
independent points uniformly drawn from [0,1]". We define a trajectory between xi 

x(0 = g[G''(^i)cos/ + G"'(;c2)sin/J (2) 
and X2 by : 

where G is the centered and reduced Gaussian distribution function. According to 
the gradual deformation algorithm described in the aforementioned French patent 
2,780,798, for any t, x(t) is a uniform point in [0,1]". When the two points are fixed, 
the t rajectory o f t he g radual m igration b etween t hem i s e ntirely d etermined. T he 
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place change of one of the two points will change the nnigration trajectory. It can be 
shown that the trajectory defined by Equation (2) is symmetrical in relation to the 
centre of the domain [0,1]", whatever the number of dimensions n. This suggests 
that, even if the two points are isolated in a corner/side of the domain, the trajectory 
of the gradual migration between them can always reach the opposite part of the 
domain. 

Migration of a non-stationarv Poisson point process 

[0035] Consider now the Poisson point process in a domain D of density X{x) 
of general shape. The number of points in D is a random Poisson variable of 
average A.{D) = j^X{x)dx . These points are independently distributed in D 

according to the probability density function : 

fix) = Mx) I Ji(D) xeD (3) 

[0036] Simulation of a Poisson point process of density X{x) in D can be 
carried out in two stages : 

generating a number n according to Poisson's law of average equal to X{D), then 

generating n points in D irrespective of one another according to the same 
probability density f(x). 

[0037] The example of Fig.lA shows a density function, and Fig.lB shows a 
realization of a Poisson point process of density X{x). 

[0038] If law f(x) is simulated by Inversion of its distribution function, then 
point X corresponds to a uniform vector u. The gradual migration algorithm can thus 
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be applied to the Poisson point process of density X{x), By construction, this 
method preserves the density and the number of points of the initial process. Figs. 
2A to 2D show four migration trajectory examples. 

[0039] Take the example of the non-stationary Poisson point process in a 
two-dimensional domain D = [0.1]^ whose density increases linearly In axis X and 
is constant in axis Y. Let (x,y) be the vector of the coordinates of a point in D of this 
process. Then, (x,y) admits the bivariable probability density : 

f(x,y) = 2x, (x,y)G[0,\f (4) 

[0040] Simulation of a point according to the above law (4) is simple : 
generating abscissa x according to the linear density f(x) = 2x, then 
generating ordinate y uniformly between 0 and 1 . 

[0041] Let : 

X = 4u 

y = v (5) 
where u and v are two independent and uniform numbers between 0 and 1. 

[0042] Thus a gradual migration algorithm (2) can be applied to the uniform 
vector (u,v) in order to establish a migration trajectory for point (x,y) in D. 

[0043] In practice, density function A,(x) often comes in form of a grid. By way 
of example, consider the case of a two-dimensional grid of M x N nodes. Let X| and 
yj be the coordinates of node (i,j). The marginal law for X| is then : 

/U) = S/(^,,y,) (6) 
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and the conditional law for yj, knowing Xj, is : 

[0044] The points of the Poisson process of density A,(x) can thus be readily 



f.^yj)-f{x,,yj)in^i) (7) 

generated on a grid from these monovariable probability laws. 

[0045] By way of illustration, a two-dimensional density function is 
constructed from a centered and reduced Gaussian simulation. The variogram is 
anisotropic and of Gaussian type. The principal anisotropy directions are diagonal 
in relation to the coordinates system. The scale factors in the principal anisotropy 
directions are 0.3 and 0.1 respectively. The size of the field is 1x1 and it is 
discretized in 1000x1000 pixels. The Gaussian numbers are changed into positive 
numbers according to the following expression : 

X{x) = 400e^^'^ (8) 

where Y(x) is the Gaussian simulation. Fig.lA shows the density function thus 
constructed. Fig. 1B shows a Poisson point process realization admitting the 
density function of Fig.lA. 2000 points are generated according to the sequential 
method. 

[0046] Figs. 3A to 3H illustrate the evolution of a Poisson process realization 
during gradual migration of the points. It can be seen that the point density is 
respected during the migration. 
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Migration of conditional Boolean simulation objects 

[0047] The above method can be immediately applied to the migration of a 
domain S of any form in D. In fact, the migration of a point in S can be carried out in 
D using the truncated probability density function : 

f,ix) = A{x)\^^,/X{S\ xeD, SczD (9) 
[0048] Using the Markovian iteration algorithm described in the reference to 

Lantuejoul, 1997, it is possible to simulate a Boolean model in a domain D knowing 

that two subsets Ci and Co of D respectively belong to the union of the objects and 

to its complement. Then, the migration algorithms in any domain can be used for 

gradual deformation of the Boolean simulations conditioned by geometric well data. 

In fact, from a conditional Boolean simulation and without compromising 

conditioning by Ci and Co, the objects must move only in their respective domains 

defined according to the shape of the objects and the configuration of Ci and Co. 

[0049] Consider an object A of a conditional Boolean realization, which 
includes a subset Cia of Ci and excludes Co- If, during migration of A, it still has to 
include Cia but excludes Co, the allowed migration domain for object A is : 

D, = {x:C,nA,=0;C,^^A,} (10) 

[0050] Fig. 4 shows a case with three conditioning points and the eight 
possible domains of migration of a disk. If, for example, the disk is supposed to 
always cover points (a) and (b) but avoid point (c), then its centre can move only in 
domain 3. 
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Migration with appearance and disappearance of points 



[0051] The number of points in D of a point process follows a Poisson's law 
of parameter X(D). It is therefore necessary to vary the number of points in D during 
their migration. In this section, two methods of constructing Poisson number chains 
are first presented, and then how to migrate between two sets of points whose 
cardinals are not identical is presented. 

[0052] A chain of Poisson numbers is attempted to be constructed between 
two Poisson numbers Ni and N2, independently generated by inversion of the 
distribution function. Let Ui and U2 be two independent uniform numbers (between 
0 and 1) from which numbers Ni and N2 are obtained. 

N2=F-\U,) (11) 

where F'^ represents the inverse distribution function of Ppisson's law. According to 
the gradual deformation algorithm, a chain of uniform numbers can be constructed 
between Ui and U2 by : 

U(0 = G[G"'(C/,)cosr + G"^(J72)sinrJ (12) 

[0053] By inversion of the distribution function of Poisson's law, a chain of 
Poisson numbers is obtained. 

Nit) = F-'[U(i)] (13) 
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[0054] The inverse distribution function can be calculated by dichotomy. 
However, if the parameter of Poisson's law is too great, this method remains quite 
costly. 

[0055] In order to save calculating the inverse distribution function of 
Poisson's law, another way of generating Poisson's numbers can be considered. It 
is well-known that the number of jumps of a Poisson process of parameter 1, in a 
segment of length X, precisely follows a Poisson's law of parameter X. Simulation of 
a Poisson process is easy to do. In fact, the intervals between two consecutive 
jumps of the process are independent and follow the exponential law of parameter 
1. By gradually deforming the exponential segments, a chain of Poisson numbers is 
obtained. It is easy to construct a chain of exponential numbers because 
calculation of the inverse distribution function of the exponential law is easy. It can 
be reminded here that the distribution function of the exponential law of parameter 
1 is written as follows : 

F(s) = l-e-\ s>0 (14) 
and its inverse function is : 

F-'(w) = -ln(l-w), 0<w<l (15) 

[0056] The fact remains that, since the number of jumps in the segment of 
length A, is A. on average, the number of chains of exponential numbers is close to X 
on average. 

[0057] Now the migration of a set of Ni points (set 1 ) to a set of N2 points 
(set 2) is considered. Since Ni ^ N2, the migration from set 1 to set 2 necessarily 
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implies the appearance or the disappearance of certain points. The number of 
points that will appear or disappear is determined by the chain of Poisson numbers 
between Ni and N2. The migration algorithm is as follows : 

a) Calculate the maximum number Nmax of chain N(t). 

b) Complete set 1 with Nmax - Ni points and set 2 with Nmax - N2 points. 

c) For each set, order the points from 1 to Nmax- All the initial points are ordered 
from the beginning. 

d) Calculate the migration trajectory of point n of set 1 to point n of set 2 

(n=1,2,....Nmax). 

e) For each state t of the set, remove the Nmax - N(t) last points. 
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